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Abstract 



ing test of models of diffractive deep-inelastic scattering at HERA. We present a comparison 
of several colour-singlet exchange models with thrust transverse momentum data from a 
recent HI analysis. We also study the effect of constraints imposed on the kinematically- 
accessible phase space by data selection cuts and find that, as a result of the pseudo-rapidity 
i-Q ' cut which is used by HI to select diffractive events, there is no dijet contribution at low 

»-H transverse momenta. We are able to fit the large transverse momentum part of the data 

a 1 with a two-gluon dijet model. The results of this analysis are compared with a previous 

Q_) | study of large pseudo-rapidity gap structure function data, and we discuss ways in which 

<"^ ■ one might reconcile the results of the two analyses. We conclude that a significant small-/3 

3-jet contribution is probably required to explain the data, and show that the combina- 
tion of a two-gluon dijet model and an exponentially-decaying thrust transverse momentum 
j^ ■ distribution provides a good fit over the entire kinematic range of the thrust data. 



1 Introduction 

With the increasing numbers of models [Q-|[ which propose to describe diffractive deep-inelastic 
scattering (DIS), it is important to find ways to test them. One method is to examine events 
with a large pseudo-rapidity gap between the proton remnant and the diffractive final state. 
The requirement of a large pseudo-rapidity gap means that the experimental cuts are selecting 
a sample with a strongly reduced phase space H . An analysis of the effects of pseudo-rapidity 
cuts in terms of diffractive structure functions has been reported previously H . A further test 
of models of diffractive DIS is to study the unintegratcd transverse momentum distribution of 
diffractive events flCl-flT 
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Here we provide a detailed analysis of the thrust transverse momentum distribution reported 
recently by HI fl3| . We study the restrictions on available phase space for diffractive inter- 
actions which are imposed by a combination of pseudo-rapidity cuts and other data selection 
cuts. We find that these cuts reject dijet (quark-antiquark diffractive final state) events with 
small transverse momenta relative to the virtual photon-proton axis. An additional contribu- 
tion is therefore required to fit the data at small transverse momenta. At the values at which 
the low-momentum dijet cut-off occurs the data also changes slope, having an approximately 
quadratic fall-off with increasing transverse momentum up to about 4-10 GeV 2 which changes to 
a quartic fall-off at larger transverse momenta. The (/3-dependent) position of the low transverse 
momentum dijet cut-off corresponds rather closely to the transverse momentum values at which 
the slope of the cross section changes. 

In this paper we consider how dijet events contribute to the transverse momentum distri- 
bution, and discuss the role that higher-multiplicity diffractive final-state events are expected 
to play. The colour-singlet exchange models of dijet production which we consider here are: 
the form factor vector pomeron models of Donnachie and Landshoff (D-L) [g, [| and Ellis and 
Ross (E-R) |j],[| (in the Feynman gauge and in a non-covariant gauge), the scalar pomeron 
model of Vermascren et al , a leading- twist single-gluon exchange model, and the two-gluon 
Landshoff-Nachtmann model described in Q. The Donnachie-Landshoff (non-covariant gauge) 
and two-gluon models provide a good fit to the quartic fall-off of the diffractive cross section 
with thrust transverse momentum, while the other models predict only a quadratic fall at large 
transverse momenta. 

The result of the thrust transverse momentum analysis appears, at first, rather in contra- 
diction with a recent analysis of large pseudo-rapidity gap structure function data ||. The 
structure function analysis appeared to strongly favour the l/pj_ models. In particular, the E-R 
and leading- twist models fit the structure function data well over a wide range of Q 2 , (3 and xp. 
The two-gluon and Donnachie-Landshoff models fall off too rapidly with decreasing (3 to fit the 
data without a substantial contribution from higher-order diffractive final states. We explore 
two possible ways to reconcile these results: the first is to consider that the models which give a 
quadratic dijet transverse momentum dependence describe the diffractive data only at large /3, 
while the small-/3 spectrum is dominated by a two-gluon-like contribution. This is a possibility 
that cannot be ruled out by present large pseudo-rapidity gap data. The other explanation for 
these results is that both sets of data might be described by the combination of a two-gluon dijet 
component and a (possibly non-perturbative) contribution which is strongly peaked at small [3 
and small thrust transverse momentum. 

The position of the low transverse momentum cut-off is largest at small values of (3. This 
means that the extra thrust transverse momentum contribution which is required to "fill in" the 
thrust transverse momentum distribution spectrum is greatest at large diffractive masses, i.e., 
at small f3. Since (3 plays the role of the scaling variable x in standard DIS, we conclude that the 
additional contribution is probably due to the gluonic structure of the pomeron. That is, the 
additional contribution needed to describe the thrust transverse momentum spectrum at small 
transverse momenta is probably due to 3-or-more parton diffractive final states, predominantly 
quark-antiquark-gluon ("3-jet") diffractive final states. Another reason for this conclusion is 
that most of the HI thrust transverse momentum distribution data corresponds to very small 
values of (3, and it is at such such values that one might expect the contribution from 3-jet final 
states to be significant @,E3, 15[. The combination of a two-gluon dijet model and a simple 
/3-dependent Gaussian fits the data well over the entire kinematic range in the HI analysis. 



The rest of this paper is organized as follows: in Sect. we present a discussion of the usual 
diffractive DIS kinematics, in which we introduce the kinematic variables which are of particular 
interest in this analysis. Sect. || contains an analysis of the effects of data selection cuts. In Sect.H 
we describe several models of diffractive dijet production and compare these models with thrust 
transverse momentum data from a recent HI study [ pjfl . In Sect, ra we compare our results with 
a recent analysis of large pseudo-rapidity gap structure function data and discuss two possible 
ways to reconcile the two analyses. This section includes a discussion of the possible contribution 
of 3-jet diffractive final states to the transverse momentum spectrum. Sect. |6J summarises the 
results presented in this paper and outlines further investigations which could be carried out to 
verify our conclusions. In the Appendix we give a brief discussion of the large pseudo-rapidity 
gap structure function analysis carried out previously ||. In this Appendix we extend the 
comparison of form factor models with large pseudo-rapidity gap structure function data to 
allow for a different form factor cut-off, A, and include fits for the form factor models calculated 
in the non-covariant gauge. 

2 Background 

2.1 Kinematics 

2.1.1 Standard Kinematics 

In the HERA electron-proton experiments 820 GeV protonsn collide with 27.5 GeV electrons or 
positrons. This corresponds to a centre-of-mass (CMS) energy of v^Tot ~ 300 GeV. In the 
HERA lab frame the positive z-axis is defined by the forward proton direction with the origin 
at the interaction vertex. We consider diffractive deep-inelastic e — P scattering, 

e( Pe ) + P(P) ^ e(p' e ) + P(P') + X(X), (1) 

where the momenta of the particles are shown in brackets. One may consider that the interaction 
proceeds by virtual photon-pomeron deep-inelastic scattering, 

1 *{q) + IP(P ] p)^X(X) 1 (2) 

where P P = P - P' . 

We use the usual kinematic variables of decp-inclastic scattering, 

Q 2 = -q 2 , X =7T5 — > andy=— - — , (3) 

2-P • q x Srot 

where Q 2 is the negative four-momentum squared of the virtual photon and x is the Bjorken 
scaling variable. We also define W 2 , the mass squared of the total hadronic system (X + outgoing 
proton), by 

W 2 = (P + q) 2 . (4) 

Additionally, for diffractive scattering we define: 
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where tp is the momentum transfer at the proton vertex and is constrained by experimental 
cuts to be small (\tp\ < lGeV 2 ), xp is the fraction of longitudinal momentum of the proton 
carried by the pomeron, and x = fix p. The mass squared of the diffractive system X is M\, 
and the proton mass is neglected. 

The pseudo-rapidity, r\, of an outgoing particle is defined in the laboratory frame in terms of 
its polar angle with respect to the proton direction by 

r/ = -lntan(-|M . (6) 

2.1.2 Exchanged Quark Virtuality, Thrust and Transverse Momentum 

For the leading-order (to order a s ) diffractive process shown in Fig.|l| we introduce a further 
invariant, the four-momentum squared of the struck quark, k 2 . In the 7* — IP CMS the virtuality 
of this quark can be expressed in terms of other invariants and the polar angle with respect to 
the 7* — IP axis, by 

fc 2 = - Q2+ Mi " (l-cosfl cms ). (7) 




Figure 1: Diffractive DIS via pomeron exchange. 

A similar expression can be formed for interactions where more than two diffractive final-state 
partons are produced. 

There are two further quantities of interest here. The first is thrust, T, which describes the 
principal axis of jet momentum and the spread of final-state particles. More formally, the thrust 
axis is defined to be parallel to the vector n which maximises the weighted sum 

1 N 

T=(— 7r — )-max^kVn|, (8) 

22i=i\Pi\ n <=i 



where the sum is taken over the whole diffractive final state and T is the thrust variableR 
The thrust axis in diffractive dijet production is the axis defined by the back-to-back final-state 
quarks in the virtual photon-pomeron CMS and is approximated experimentally by the thrust 
axis described by the (hadronic) diffractive final state. For a dijet event T is equal to 1. while for 
a symmetric 3-jet event T = 2/3 and a completely isotropic multi-jet event will have T = 1/2. 
In general, T can take any value between 1/2 and 1. 

Another useful kinematic quantity is the thrust transverse momentum, p\, defined in the 
virtual photon-pomeron CMS. For events with dijet final states, p± is simply the component 
of momentum of one of the final state hadronic jets which is transverse to photon-proton axis 
(equivalently, to the photon-pomeron axis[j). This is shown in Fig.g. In the more general case of 
final states with two or more partons, the experimentally measured transverse momentum vari- 
able is the transverse momentum of the thrust axis with respect to the photon-proton direction. 
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Figure 2: Definition of transverse momentum variable p\ (in 7* — IP CMS). 



The transverse momentum variable defined above is different to that defined in |l6[ where 
phase space constraints due to pseudo-rapidity cuts |l|,|9[] are discussed in terms of a lower limit 
on the transverse momentum of the final-state parton which couples to the pomeron. The effect 
of pseudo-rapidity cuts on the 3-jet thrust transverse momentum distribution is more subtle as 
there exist 3-jet final-state configurations in which Pj_ thrust = 0, for example, but for which the 
events will survive most pseudo-rapidity cuts. Events where there is a parton emitted far in 
the forward direction, however, are expected to be rejected by pseudo-rapidity cuts and hence 
there is likely to be a systematic reduction in the low-pj_ 3-jet thrust transverse momentum 
spectrum. Contraints on thrust transverse momentum due to data selection cuts could, at best, 
be investigated through Monte Carlo studies of diffractive DIS. 

3 Phase Space Analysis of Data Selection Cuts 

We start this section with a summary of the data under investigation and then discuss the effects 
of the cuts used to select the diffractive DIS sample. 

HI measured the diffractive DIS cross section as a function of thrust transverse momentum, 
Pj_, and final-state diffractive mass, Mx- They presented their results in the form [ O 

2 Some references use a slightly different definition of T. The normalisation of T we have adopted here is that 
used in [113]. 

3 Here we assume that, for small tp, the pomeron is emitted parallel to the proton direction. 



I dV ( 2, 



(9) 

fixed M x 



The diffractive cross section has been integrated over Q 2 and xp in the range 



2 



Q z : 10 -> 100 GeV z , 

xjp : 10~ 4 -> 3.2 x 10~ 2 , (10) 

and the data was binned with average diffractive masses of 

M x = 7.01, 9.40, 12.82, 16.85, 21.20, 28.58 GcV. (11) 

This corresponds to integrating over small values of /3, for example: for Mx = 7.01 GeV, over 
the range of Q 2 in this analysis (3 ranges from 0.17 to 0.67; whereas for Mx = 28.58 GeV, (3 
ranges from 0.01 to 0.11. Here "small" is to be interpreted as meaning values of j3 for which one 
might expect 3-jet processes to provide a significant contribution to the diffractive final state 
(specifically, /3 < 0.3 §). 

The HI transverse momentum distribution data p| has a l/p\ fall-off at low transverse 
momenta, steepening to a l/p\ fall-off at higher transverse momenta. This change in slope 
occurs between about 1-10 GeV 2 for the values of Mx in this study. The l/p\_ tail in the data 
is not a phase space effect, as this tail is present for all values of Mx a long way away from the 
phase space limit oip 2 lmax = M x /A. For Mx — 28.58 GcV, for example, the data changes slope 
at about p\ = 10 GeV 2 . 

In an analysis of this data one must take account of the effects of experimental cuts on the 
space of Q 2 and xp which was integrated over. There are two cuts which restrict the accessible 
range of Q 2 and xp for a given final-state transverse momentum and diffractive mass: 

1. Cut on y: y < 0.5. 

This cut is imposed to restrict the photo-production background, as well as removing events 
where the diffractive system is strongly boosted in the backward direction [^3|. Using the 
relation (0) between y and the other kinematic variables we get the restriction on xp, 

2{Q 2 + M 2 X ) 

xp > x , (12) 

J Tot 

which corresponds to a lower limit on the integral over xp as a function of Q 2 at fixed Mx- 

2. Pseudo-rapidity cut: require pseudo-rapidity gap greater than ?7 max = 3.2. 

It was shown in P,[L6| that imposing a pseudo-rapidity cut removes part of the low-pj_ dijet 
contribution to diffractive DIS. For a given range of kinematic variables Q 2 , M x and xp, a 
pseudo-rapidity cut corresponds to imposing a lower transverse momentum cut-off, p 2 j_ min , on 
accepted events. Therefore, when calculating the cross section at a given p\, say Pj_ fixed , one 
should only integrate over the region of Q 2 and xp for which p 2 Lmin < plfi xod - The region of 
phase space for which pl^^ > Pj^fixod wu l have been removed by the pseudo-rapidity cut. 
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Figure 3: Phase space plot for p^_ min = 4, n niax = 3.2 (hadron-level cut), and Mx = 16. The 
horizontal lines show upper and lower limits of xp in data sample, while the upper slanted line 
shows the upper limit on xp as a result of pseudo-rapidity cuts, and the slanted lower line shows 
the lower limit of xp as a result of the y cut. 



The above discussion of the effect of pseudo-rapidity cuts on available phase space is illus- 
trated in Fig.H. In this figure, the region in the Q 2 — xp plane over which one would integrate 
to extract the cross section at p\ = 4 GeV 2 is shown for a diffractive mass of Mx = 16 GeV and 
a pseudo-rapidity cut of 3.2 (typical HERA values). The upper shaded area shows the region 
excluded by the pseudo-rapidity cut. In this region of Q 2 and xp dijet events with p\ < 4 GeV 2 
are removed by the pseudo-rapidity cut. The lower shaded area shows the region excluded by 
the cut on y. As a result, there is a restricted region of phase space populated by dijet events 
with thrust transverse momentum of 4 GeV . 

For large diffractive masses and small p\, the data selection cuts exclude the entire Q 2 — xp 
plane, leading to a sharp low-momentum cut-off in p\ below which there are no values of the 
kinematic parameters Q 2 and xp in the HI study for which dijet events can contribute to 
the diffractive cross section. The position of the cut-off depends on Mx, and moves to larger 
transverse momentum values for larger diffractive masses. In Fig.H, the position of the low- 
momentum dijet cut-off is shown for three diffractive masses: Mx = 28.58 GeV, 21.20 GeV, and 
16.85 GeV. For smaller diffractive masses the transverse momentum distribution of dijet events 
is not strongly restricted by the data selection cuts and, for Mx < 15 GeV, there is no lower p± 
cut-off from the data selection cuts. 

This analysis of the effect of data selection cuts on the dijet contribution to diffractive DIS 
therefore demonstrates that an additional contribution is needed to describe the HI thrust trans- 
verse momentum distribution data, particularly at low transverse momenta. 



4 Analysis of Diffractive Dijet Production Models 

In this section we present the fits of various dijet models to the HI thrust transverse momentum 
distribution data, starting with a brief description of the diffractive dijet models under consid- 
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Figure 4: Position of dijet transverse momentum cut-offs at different diffractive masses (two- 
gluon model). The solid line corresponds to a diffractive mass of 28.58 GeV, while the dot-dashed 
line is for 21.20 GeV, and the dotted line represents a diffractive mass of 16.85 GeV. 

eration. The contribution of higher-order diffractive processes, such as quark-antiquark-gluon 
diffractive final states, is discussed in the following section. 

4.1 Models of Diffractive Dijet Production 

The first type of model we studied is the form factor approach of Donnachie and Landshoff [0,p[ 
and Ellis and Ross 0, H] in which the diffractive interaction is assumed to proceed by the 
exchange of a C— even, vector-like pomeron. This approach allows for a direct coupling between 
the pomeron and off-shell partons. In order to maintain the experimentally-observed scaling of 
the diffractive structure function, a form factor is inserted to soften the virtual quark-pomeron 
vertex. In the original phenomenological model of Donnachie and Landshoff 0,0], applied to 
deep-inelastic scattering, one inserts a form factor at the virtual quark-pomeron vertex of 



f(k 2 ) 



A 2 



A 2 



(13) 



Here k 2 is the four-momentum squared of the exchanged quark coupling to the pomeron (the 
other quark at this vertex is near mass shell and hence its momentum may be neglected), and 
A 2 is of the order of the hadronic scale (A 2 ~ 1 GeV 2 ). This serves to cut off contributions at 
large parton virtualities. In the modification suggested by Ellis and Ross |j|], one chooses the 
form factor 



!{k 2 ) 



A 2 



A 2 - fc 2 ' 



(14) 



This choice of form factor leads to a cross section with a logarithmic dependence on k 2 , which 
therefore receives contributions from the entire k 2 spectrum]]. This model is motivated by 

4 Both form factor models reduce to the original Donnachic-Landshoff model of hadron-hadron scattering 
]17J-|20| for the exchange of a pomeron between two on-shell partons. 



the observation that diffractive DIS is relatively insensitive to pseudo-rapidity pi|-j23|], as well 
as the correlation between between exchanged quark virtuality and pseudo-rapidity which was 
demonstrated in 0|. 

The other models we consider were described in detail in H. The first is a model in which 
the "pomeron" is described by the colour-singlet part of single-gluon exchange, that is, the 
exchange of one hard gluon with subsequent exchanges of soft gluons which "dress" the gluon 
colour enabling the formation of a pseudo-rapidity gap. We also consider a direct-coupling scalar 
pomeron model j|] . Finally, we consider the two-gluon Landshoff-Nachtmann model p] following 
the analysis of Dichl || . In this approach the colour-singlet state is modeled by the exchange of 
two non-interacting gluons which are described by non-perturbative propagators and which are 
allowed to couple to the two quarks in all possible combinations. 

It should be noted that the form factor models are not gauge-invariant @,0]. The sum 
of the diagrams which contribute to two-gluon exchange, as shown in Fig.|5j, is an explicitly 
gauge-invariant quantity. Comparing the two approaches, the diagrams where the colour-singlet 
exchange is modeled by a hadron-like "pomeron" exclude the two diagrams in which the two 
gluons do not couple to the same quark line. At best, therefore, one might hope for the form 
factor models to be physically meaningful in a gauge where the contribution from these two 
diagrams is minimal. Diehl [J7| has suggested a non-covariant gauge in which the D-L model 
mimics many of the features of the manifestly gauge-invariant Landshoff-Nachtmann two-gluon 
model. We find that studying the unintegrated transverse momentum distribution highlights the 
gauge non-invariance problem with the form factor approach, and we have chosen to calculate 
the form factor models in the non-covariant gauge in this analysis. We have also repeated the 
structure function analysis of M for the form factor models in this gauge (see Appendix). 





Figure 5: Two-gluon exchange graphs used to model pomeron exchange. 



4.2 Transverse Momentum Distribution Fits with Dijet Models 

The diffractive scattering cross section is calculated by integrating the cross section from the 
colour-singlet exchange models over the region of phase space in Q 2 and xp accessed in the HI 
analysis. The common feature of all the fits is that, as a result of the large pseudo-rapidity gap 
requirement and the cut on y, the dijet contribution is cut off sharply at low p^. For example, for 
a diffractive mass of Mx = 28.58 GcV, the phase space restrictions affect the dijet contribution 
below about 10 GeV 2 , completely removing the contribution from events with p\ < 4 GeV . For 
Mx = 16.85 GeV, the turn-over is much sharper and the spectrum cuts off between 1-2 GeV 2 . 



The curves in Figs. [6] and [?] represent the dijet model fits to the HI data by fitting only the 
part of the curves after the phase space-induced turn-over. This shows how the models describe 
the large-p^ behaviour of the cross section. In each graph the HI data [Q is shown as solid 
points with error bars representing statistical and systematic errors added in quadrature. 

The fits shown in Figs.|6(a) and 7(a) arc of the Donnachie-Landshoff form factor model in 



the Feynman gauge (solid line) and in the non-covariant gauge (dashed line). In both cases we 
have assumed a pomeron intercept of 1.2 and a low-momentum cut-off in the form factor of 
A = 1.2 GeV 2 . These graphs show that the Donnachie-Landshoff model in the non-covariant 
gauge provides a good fit to the large-pj_ part of the transverse momentum spectrum, that is, 
it predicts a \/p\_ fall-off for p\ > 10 GeV 2 . The Feynman gauge curve gives a poor fit, mainly 
due to a large-pj_ tail w hich is not present in the non-covariant gauge approach M . 

Figs. 6(b) and |7(b)| show fits of the Ellis-Ross form factor model to the transverse momen- 



tum distribution data with ap(0) — 1.2 and A = 1.2 GeV 2 . Again, a high-p^ tail is evident 
in the cross section calculated in the Feynman gauge. The E-R model fails to describe the 
large transverse momentum behaviour of the data in either gauge, predicting a l/p\ fall-off at 
large transverse momenta. The scalar pomeron and single-gluon exchange models are shown in 
Figs. 6(c) and [7(c) . These models also fail to describe the pj 4 dependence of the data at large 



p\ , giving only a p , 2 fa ll- off. 

Finally, in Figs. |6(d) and |7(d)| we study the two-gluon model of Diehl ||. This model fits 



the data well above the phase space-induced low-p^ cut-off. The fit here is similar to the fit 
obtained using the Donnachie-Landshoff form factor model. 

5 Comparison of Structure Function and Thrust Trans- 
verse Momentum Distribution Analyses 

The results of fitting the colour-singlet exchange models to the thrust transverse momentum 
distribution data are surprising when compared with the results from the comparison of these 
models with the large pseudo-rapidity gap structure function data in [p| . The structure function 
data is obtained from the p^-integrated cross section. The Ellis-Ross and singlc-gluon exchange 
models provided an acceptable fit to the large pseudo-rapidity gap structure function data, but 
did were unable to describe the shape of the thrust transverse momentum distribution data. By 
contrast, the Donnachie-Landshoff and two-gluon models fit the thrust transverse momentum 
distribution data well, but under-estimated the structure function data badly at low /3. There 
are two possible ways to reconcile these results. The first possibility is that E-R and single-gluon 
exchange models may describe the data at large f3 (in the region where the operator product 
expansion is believed to be reliable in describing diffractive DIS), but that the leading-twist 
description fails at small (3. The other possibility is that the D-L and two-gluon models are able 
to fit both sets of data with a substantial contribution from multi-jet diffractive final states. This 
latter possibility requires the multi-jet contribution to the low-momentum p\ data to be most 
significant at small-/? and to have a thrust transverse momentum distribution which is strongly 
peaked near zero. 

We discuss these possibilities in more detail in this section. 



10 



a dp 2 




20 40 60 80 100 



pi GeV 2 



a dp 2 ± 

♦ 

0.1 t 




20 40 60 80 100 



v\ GeV 2 



(a) Donnachie-Landshoff form factor model in 
non-covariant gauge (solid line) and in Feyn- 
man gauge (solid line). Here A = 0.2 GeV 2 
and a P (0) = 1.08. 



(b) Ellis-Ross form factor model in non- 
covariant gauge (solid line) and Feynman 
gauge (dashed line). Here A = 0.2 GeV 2 and 
a P {0) = 1.08. 
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(c) Single-gluon exchange model (solid line) 
and Scalar Pomeron exchange model (dashed 
line). 



(d) Two-gluon exchange model. 



Figure 6: Fits of colour- singlet exchange models to transverse momentum distribution data from 
M. The fit is made to points after the turn-over of the curves (see text for explanation) . The 
curves correspond to the models, while the solid points are HI data. The plots show the fits for 
diffractive masses in the region of Mx = 28.58 GeV. The statistical and systematic errors have 
been added in quadrature. 
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(a) Donnachie-Landshoff form factor 
model in non-covariant gauge (solid line) 
and in Feynman gauge (dashed line). Here 
A = 0.2 GeV 2 and a P (0) = 1.08. 



(b) Ellis-Ross form factor model in non- 
covariant gauge (solid line) and in Feyn- 
man gauge (dashed line). Here A = 
0.2 GeV 2 and ct P (0) = 1.08. 
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(c) Single-gluon exchange model 

line) and Scalar Pomeron exchange model 

(dashed line). 



(d) Two-gluon exchange model. 



Figure 7: Fits of colour- singlet exchange models to transverse momentum, distribution data from 
\19jj. The fit is made to points after the turn-over of the curves. The curves correspond to the 
models, while the solid points are HI data. The plots show the fits for diffractive masses in the 
region of M x = 16.85 GeV. 
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5.1 /5-Dependent Structure Function Models 

It is possible that the Ellis- Ross or leading-twist single-gluon exchange models describe a physical 
picture of the data for large /?, but that some other (possibly higher-twist) component is required 
at low /?. The E-R model provided a good fit (x 2 /dof = 37/41) to the large pseudo-rapidity gap 
structure function data studied in H and described the data well over a wide range of Q 2 , j3 and 
xp . Most of the large pseudo-rapidity gap structure function data is at low /?, and the transverse 
momentum distribution data is almost all at low /?, however, and there is not sufficient large-/? 
data available to test this possibility. 

The present thrust transverse momentum distribution analysis demonstrates a need for a 
substantial extra non-dijet contribution to the cross section. Due to the low-pj_ dijet cut-off. 
some additional contribution is needed to explain the low transverse momentum data. The 
position of the cut-off shows also that the size of this extra contribution is /3-dependcnt, and 
peaked at low /?. If this low-/? contribution is included in the E-R model fits to large pseudo- 
rapidity gap structure function data it will cause the E-R model to overshoot the structure 
function data at low /?. The need for a substantial low-pj_, low-/? contribution to explain the 
thrust transverse momentum distribution data therefore presents a major challenge to E-R and 
leading-twist models. 

5.2 Two-Gluon Dijet Contribution Plus 3-Jet Contribution 

The M^-dependent position of the dijet low transverse momentum cut-off is such that at large 
diffractive masses the region where dijet events cannot contribute to the diffractive scattering 
cross section is larger than for smaller diffractive masses. This dominance at small /? is what 
is expected from quark-antiquark-gluon diffractive final states &H 111]- It seems reasonable 
therefore, given also the very small values of /? in the HI thrust study, to conclude that the 
extra contribution required to explain the thrust transverse momentum distribution spectrum is 
probably provided by 3-jet diffractive final-state events. This conclusion is consistent with |8|, 
where it was estimated that the 3-jet contribution dominates the diffractive structure function 
for (3 < 0.2, and provides a significant contribution to F 2 for /? < 0.4. This hypothesis is 
also consistent with the HI findings that the average thrust of diffractive DIS events is less than 
0.9 g|. 

If the 3-jet contribution is indeed dominant at small (3 (large Mx), one would expect to 
see the pattern of diffractive structure function results found in || for the D-L and two-gluon 
models. In this case the dijet contribution is expected to under-estimate the diffractive structure 
function at small (3. 

Here we discuss the low-/? component needed to give a good fit to the thrust transverse 
momentum distribution data. We assume in this discussion that the additional contribution 
at low-/? originates from quark-antiquark-gluon diffractive final states. Considering the results 
shown in Figs. |6| and [[[ and the analysis in |9j], one can immediately make two hypotheses about 
the distribution of "3-jet" processes in the diffractive sample: 

• 3-jet events dominate the large pseudo-rapidity gap sample at low (3: 

The requirement of a large pseudo-rapidity gap between the proton direction and the diffrac- 
tive final state removes the low-p"^ dijet events from the sample ||[16|], hence higher-multiplicity 
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final states are required to explain the low transverse momentum part of the spectrum. Fur- 
thermore, comparing Figs.o and ffl, we see also that the cut-off in the dijet contribution from 
the phase space restrictions is shifted to smaller p\ for smaller values of Mx- This is to be 
expected as smaller Mx corresponds to large values of /3, for which there is a weaker restriction 
on the dijet phase space &pl|. Also the ratio of multi-jet events to dijet events is expected to 
be smaller at large f3 || . 

• The thrust transverse momentum distribution of 3-jet events is strongly peaked about 0: 

This conclusion is partly based on the success of the two-gluon dijet model in fitting the pj 
behaviour of the cross section at large transverse momenta, and also on the inability of all the 
dijet models to describe the low-p^ part of the spectrum due to data selection cuts. If the 3-jet 
events were to describe the part of the transverse momentum spectrum which is not described by 
dijet events, one would also require them to have thrust transverse momentum strongly peaked 
about p\ = 0, in order to not effect the spectrum at larger transverse momenta. 

A strongly peaked thrust transverse momentum distribution might be expected from a two- 
gluon "partonic" model of the pomeron. In such a model, the leading-order process is one where 
the photon interacts with one of the gluons in the pomeron, producing a diffractive quark- 
antiquark pair. The third final-state particle is the other constituent gluon in the pomeron, 
the pomeron remnant, which would be traveling in approximately the same direction as the 
initial proton. Therefore, especially if there is a reasonably isotropic distribution of momentum 
between the two constituent gluons in the pomeron, one would expect the thrust axis of the final 
state to be rather closely aligned with the proton-photon axis. This hypothesis does, however, 
require a /3-dependent non-perturbativc spreading of the transverse momentum distribution of 
the remnant gluon of up to Pj_= 6 GeV . The intrinsic transverse momentum of the pomeron has 
been discussed recently in pi| . It is also possible that the peaked thrust transverse momentum 
distribution might be due to the non-zero momentum transfer at the proton vertex [E5| . It is not 
possible to explore this latter possibility further at present as there is insufficient data describing 
the ijp-dependence of diffractive DIS. 

5.2.1 Parameterization of the Extra Thrust Contribution 

We fit the thrust transverse momentum distribution data with a combination of the two-gluon 
model and a parameterization for the "3-jet" contribution to diffractive DIS in which the invariant 
amplitude of the thrust transverse momentum distribution is described by the GaussianH: 



/( ° 2 ^ = o¥ e — (15) 

This provides a /3-dependent distribution which falls sharply with increasing p\. Since the 
HI data contains no information about the absolute normalisation, in these fits the overall 
normalisation of each contribution is a free parameter, giving two adjustable parameters when 
combining the two contributions. The most unsatisfactory feature of the "3-jet" interpretation of 

5 Note that since the definition of the transverse momentum variable in this analysis is different to that used in 
the diffractive structure function study in |H|, the parameterization used to fit the thrust transverse momentum 
distributions here cannot be directly applied to the structure function analysis. 
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this extra contribution is the need for an exponential /^-dependent cut-off in the thrust transverse 
momentum distribution. 

The combination fit is shown in Fig.H. This fit describes the transverse momentum dis- 
tribution data well in the range of diffractive masses studied by HI, only underestimating the 
data slightly at Mx = 16.85 GeV. This confirms that the data may be described by the com- 
bination of a dijet model with a quadratic fall-off with transverse momentum and a "multi-jet" 
parameterization with an exponentially decaying thrust transverse momentum distribution. 

A more complete test of this hypothesis would be to study data for which the overall normal- 
isation of the cross section at each diffractive mass is measured, as this would enable allow the 
relative contributions of 2- and 3-jet final states to the diffractive cross section to be measured. 

6 Discussion and Conclusions 

We have studied the transverse momentum distribution of diffractive DIS and compared the 
dijet component of several diffractive DIS models with recent HI data [Dj ]. 

The general trend in the data is towards a l/pj_ fall-off at low transverse momenta, steepening 
to l/pj_ at about 4-10 GeV 2 . The position of the change in slope depends on the diffractive mass, 
with the change in slope occurring at larger values of the transverse momentum variable for larger 
M x . 

The data selection cuts, particularly the pseudo-rapidity cut imposed by HI to select diffrac- 
tive events, lead to a restriction on the available phase space for a wide of the parameter space 
covered by HERA. The main result is that, for most diffractive masses and for the range of 
Q 2 and xp in this analysis, data selection cuts are expected to remove all dijet events with 
small transverse momenta. As a result, one requires a further contribution, for example quark- 
antiquark-gluon final states, to explain the small-p^ region of the HI transverse momentum 
distribution data. The region in p± where the data changes slope corresponds rather closely to 
the region where phase space effects due to data selection cuts remove the \ow-p± dijet contri- 
bution. 

We compared four classes of colour-singlet exchange models with the HI data: two form 
factor vector pomeron models, a leading-twist single-gluon exchange model, a scalar pomeron 
model, and a two-gluon exchange model. The cross section from the form factor models was 
calculated in both the Feynman gauge and in a non-covariant gauge. The difference between 
the two gauge choices was particularly evident in the large-p^ tail which one obtains using the 
Feynman gauge, which is not seen in any other approach, and which fails to describe the present 
data. The Donnachie-Landshoff form factor model (in non-covariant gauge) provided a good fit 
to the large transverse momentum part of the spectrum for all values of diffractive mass beyond 
the point where the slope changed to a l/pj_ fall. The Ellis-Ross form factor model, as well as 
the single-gluon and scalar pomeron models, predicts a l/p\ fall at large transverse momenta 
which is not consistent with the HI data. The two-gluon model fit the large-pj_ part of the 
transverse momentum distribution spectrum well. 

Since the region where dijet events cannot contribute to the diffractive cross section is much 
larger for larger diffractive masses (i.e. smaller j3), it was concluded that the extra contribution 
required to describe the HI transverse momentum distribution data was probably from 3-jet 
diffractive final states. It was found that a combination of the two-gluon dijet model and a 
simple parameterization for the extra contribution using /^-dependent Gaussian gave a good fit 
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Figure 8: Comparison of the combination of a two-gluon dijet model and a 3-jet parameterization 
with thrust transverse momentum distribution data from US] at different diffractive masses. 
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to the entire thrust transverse momentum distribution spectrum. 

The 3-jet interpretation of this extra contribution is consistent with the earlier structure 
function analysis of large pseudo-rapidity gap diffractivc DIS |9J . In this interpretation a large 
3-jet diffractive contribution is required at small j3 in addition to the dijct contribution from 
the Donnachie-Landshoff or two-gluon models. In addition, however, the present analysis also 
requires the 3-jet diffractive events to have thrust transverse momentum strongly peaked about 
zero. This analysis also suggests that the apparent success of the Ellis-Ross, singlc-gluon and 
scalar pomeron models in describing the diffractive structure function may be due to these models 
over-estimating the large-p^_ dijet contribution. 

An alternative possibility is that the E-R or leading-twist single-gluon models may provide 
a reasonable fit to the diffractive DIS cross section at moderate to large values of (3. Present 
pseudo-rapidity gap thrust transverse momentum distribution and structure function data in 
which the requirement of a large pseudo-rapidity gap places strong restrictions on the available 
phase space corresponds mostly to small /3 and at present there is insufficient data, particularly 
at large-/3, to either prove or disprove these hypotheses. 

This analysis demonstrates the worth of studying the transverse momentum dependence in 
distinguishing between models. Clearly it is of importance to obtain an improved data set in 
the diffractive scattering regime. In particular, it would be helpful if smaller xp data were 
available as the present jet production data has relatively large xp where the purely diffractive 
interpretation may be invalid. 

Further investigation is required to test the hypotheses presented in this paper. On the 
theory side, there is a need for a full 3-jet diffractive DIS calculation or at least, as far as testing 
this analysis is concerned, for a proper Monte Carlo treatment of thrust transverse momentum 
distribution of 3-jet diffractive events with proper account taken of experimental data selection 
cuts. Better statistics and finer binning in the thrust transverse momentum variable would allow 
one to explore the interface between diffractive dijet production and the low-p^ contribution. 
In a similar vein, analysing data corresponding to larger values of /?, where phase space effects 
due to pseudo-rapidity cuts are less significant and dijet events are expected to dominate over 
higher-multiplicity diffractive final states, one would expect to see the quartic p^-dependence in 
the present HI data extend to smaller values of the p±. On the other hand, restricting to smallcr 
values of /3, one would be able to push the position of the change to a quartic fall-off to larger 
values of p\ . 

A more interesting experiment would be to keep the same range of kinematic parameters 
in Q 2 , P and xp and impose a strongly pseudo-rapidity cut. This would push the position of 
the low-momentum dijet cut-off to larger values of p\. The only effect on the sharply-peaked 
low thrust-pj^ contribution would be a decrease in the overall normalisation. With fine enough 
binning in p± this should therefore lead to the production of a "dip" in the cross section in the 
p\ region between where the exponential low-p^ contribution cuts off and where the l/p^_ dijet 
contribution starts. 

Information about the overall normalisation of the cross section would enable the relative 
contributions of dijct and 3-jet events to diffractive DIS to be measured. 
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Appendix. Structure Function Analysis of Form Factor 
Models in Non-Covariant Gauge 

This paper extends the form factor analysis of |H| to the non-covariant gauge choice described 
in M . It is interesting to see whether the form factor models describe the large pseudo-rapidity 
gap structure function data in this gauge. We have also repeated the structure function analysis 
for the form factor models with a low-momentum cut-off of A = 1.2 GcV 2 in the form factor, as 
suggested by fits to vector meson production data by Donnachie and Landshoff H, rather than 
the smaller value of A = 0.2 GeV 2 used in the original pseudo-rapidity gap structure function 
analysis 0] . We compared the Donnachie-Landshoff and Ellis- Ross form factor models with the 
large virtuality constraint data from |2q] . Fitting to the 42 points considered in M , we obtained 
the x 2 parameters shown in Tablefi]. 





X 2 /dof 


a P {0) = 1.08 a P (0) = 1.2 




39/41 (37/41) 57/41 (54/41) 

81/41 (102/41) 79/41 (97/41) 
41/41 (41/41) 56/41 (55/41) 
69/41 (86/41) 68/41 (82/41) 


E-R: f(k-) ^^ fc2 

D-L: f(k 2 ) - A t k 2 
E-R: (non-covariant gauge) 
D-L: (non-covariant gauge) 



Table 1: Result of the fit of form factor models to large pseudo-rapidity gap diffractive structure 



function data from ]2w. The fits assume a cut-off in the form factors of A = 1.2 GcV . Results 



of fits with a form factor cut-off of A 2 
analysis in M are shown in brackets. 



0.2 GeV in the form factors, as was used in the original 



We can see from Table pi that changing the low- momentum form factor cut-off, A, and calcu- 
lating the cross section in the non-covariant gauge make no real difference in the fits obtained in 
the Ellis-Ross model. The most significant effect on the Donnachie-Landshoff fits is from using 
the larger form factor cut-off. The D-L model, calculated in the non-covariant gauge with a 
form factor cut-off of A 2 = 1.2 GeV 2 , gives a moderately successful fit to the diffractive structure 
function data. The principal reason that the results for form factor cut-offs A 2 = 0.2 GeV or 
A 2 = 1.2 GeV are similar is that much of the data we are studying corresponds to an exchanged 
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quark virtuality of k 2 > a few GeV 2 in the form factors, hence the low-momentum cut-off in 
the form factor has little effect. The difference in the cross sections between using the Feyn- 
man gauge or the non-covariant gauge is visible when the unintegrated transverse momentum 
distributions are studied, as can be seen in the graphs in Sect.H. 
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